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ABSTRACT 

A sample of very high resolution cosmological disk galaxy simulations is used to investigate the 
evolution of galaxy disk sizes back to redshift 1 within the ACDM cosmology. Artificial images in 
the rest frame B band are generated, allowing for a measurement of disk scale lengths using surface 
brightness profiles as observations would, and avoiding any assumption that light must follow mass as 
previous models have assumed. We demonstrate that these simulated disks are an excellent match to 
the observed magnitude - size relation for both local disks, and for disks at z=l in the magnitude/mass 
range of overlap. We disentangle the evolution seen in the population as a whole from the evolution 
of individual disk galaxies. In agreement with observations, our simulated disks undergo roughly 1.5 
magnitudes/arcsec 2 of surface brightness dimming since z=l. We find evidence that evolution in the 
magnitude - size plane varies by mass, such that galaxies with M* > 10 9 M© undergo more evolution 
in size than luminosity, while dwarf galaxies tend to evolve potentially more in luminosity. The disks 
grow in such a way as to stay on roughly the same stellar mass - size relation with time. Finally, due 
to an evolving stellar mass - SFR relation, a galaxy at a given stellar mass (or size) at z=l will reside 
in a more massive halo and have a higher SFR, and thus a higher luminosity, than a counterpart of 
the same stellar mass at z=Q. 

Subject headings: galaxies: evolution — galaxies: formation — methods: N-Body simulations 



1. INTRODUCTION 



Recent 
GOODS, 



observational surveys (e.g., COSMOS , 



Giavahsco et al.l l2004t iRix et al.l 120041 : iKennicutt et all 
20031: lYork et al.ll2000l ) have allowed us for the first time 
to statistically explore issues of galaxy formation. One 
of the immediate challenges to disk formation theory 
that is presented by these observations comes in the 
form of evidence that there has been little change in the 
sizes of disk galaxies since z—1, despite an expectation 
that disks should be growing in size over this time. 
Observations tell us that there exist large disks by z—1 
(before the universe was even half of its present age), 
and suggesting that these disks must be assembled prior 
to this epoch (e.g.. iVogt et al.l 119961: iRoche et al.l 119981: 
Lilly et ali 1998; Simar d et al.l[ l999; Ravindranat h et aLl 
20041 iFerguson et al.1 120041: iTruiillofe Aguerril 120041: 
Barden et al.ll2005t iSargent et al.ll2007l : iMelbourne et aD 
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120071 iKanwar etHI 120081 ). Recent observations 
lend strong support to the e xistence of disks 
back to even higher redsh ifts (ILabbe et al.l 
Forster Schreiber et al.l 12006 : iShapiro et al ' 



2003 



2008 



Stark et al.ll2008t I Wright et al.ll2009t I Jones et al.ll200ll 

and for disk assembly at even higher redshifts in the 
form of clump-chain galaxies in th e Hubble Ultra Deep 
Field (|Elmegreen et al.ll2(M l2007h . 

Results from large surveys suggest that there has been 
only weak evolution in disk sizes since z—1. There is 
no evidence that the size function of disks has evolved 
back to z=l dLillv et al.lll998tlRavindranath et al.ll200l 
IKanwar et al J 120081 ) . IKanwar et al.1 (|2008f ) found that 
the shape of the size function did not evolve with red- 
shift, though the normalization (or amplitude) did. This 
can be interpreted in two ways. First, the normaliza- 
tion will vary if the number density of disks varies with 
time. However, observations suggest that the number 
density of disks is constant to z—1 (jLillv et al.l 119981 : 
ISargent et al.l l2007). Second, if it is assumed that galax- 
ies dim with time (e.g., due to a passively evolving stellar 
population or decline in star formation rate), a 1 to 1.5 
magnitude dimming since 2=1 could explain the change 
in normalizatio n, while keeping the size distribution of 
disks constant ()Kanwar et al.ll2008D . Hence, there is no 
immediate evidence for a change in the sizes of galaxy 
disks with time. 

Studies of the evolution in the magnitude - size rela- 
tionship for disk galaxies back to z—1 have found sim- 
ilar results. Importantly, these studies must use sur- 
face brightness evolution to interpret evolution within 
the magnitude - size plane. This is not straightforward, 
as an increase in disk scale lengths at a fixed magni- 
tude between z—1 and z—0 can be mimicked by a de- 
crease in luminosity at a fixed size ([Truiillo fc Aguerril 
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2004). Selection effects are difficult to disentangle, 
and studies that require high and low z samples to 
adhere to the same selection biases have found that 
large disks are consistent with no evolution in sur- 
face brightness (|Simard et al.lll999t IRavindranath et all 
|2004f ). Later studies concluded that this work was 
too restrictive, and that a careful treatment of com- 
pleteness as a function of redshi ft supports s urface 
brightness d imming over time flSchade et al.l [19961: 

P^Tno o+ al I ITOOSl T ?llir of al I 11 QQBl- iRmiiiranc Hr Q^U 



Roche et all IT998I: Hillv et al I I1991 iBouwens fe Silkl 
2002t IRavindranath et al.l I2004I: ITrujillo fe Pohlenll2005l: 
Barden et al.ll2005HMelbourne et al.ll2007l: IKanwar et al.l 



2008D . Most studies conclude that pure size evolu- 
tion at a fixed magnitude is ruled out (|Melbourne et al.l 
2007). Instead, the observed surface brightness evo- 
lution is best explained by luminosity dimming, with 
the amount of dimming dependent on galaxy size so 
that lower mass galaxies have un dergone more dimming 
since z=l than mas sive galaxies ((Melbourne et al.ll2007t 
IKanwar et al.ll2008() . However, some amount of size evo- 
lution can't be ruled out, and at least weak size evolution 
is favored (IReshetnikoy et all I2003t IBarden et all I2005t 
ITrujillo fe Pohlenll2005f )~ 

While magnitude is assumed to scale with the mass 
of a galaxy, it is obvious from the above discussion that 
magnitude at a given mass is not necessarily constant 
in time. The more fundamental relation is between 
size and stellar mass. Magnitude at a given stellar 
mass is likely to increase back to z—1 due to the fact 
that there is an evolving star formation rate (SFR) - 
mass relation dGavazzi fc Scodeggiolll996t iBoselli et all 



2001; Kauff mann et all 120031: iBrinchmann et al 



Feulner et al.l 120051: l Erb et al.l 120061: ISalim et a! 



2004; 



2007 



Noeske et al.l l2007bllal: lElbaz et all 12007 : Daddi et all 
Schiminovich et all 120071: iCowie fc B arger 
Damen et al.l l2009bUal: 
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Pannella et all 12 009: _ 
Dunne et al.l 120091: iRqdighiero et al.l 120101: lOliver et al.l 



20101: iMannucci et all 120101: lLara-Lopez et al.l I2010D . 
Barden et al. I ([2005D examined the stellar mass - size 
relation for disk dominated galaxies back to 2=1, 
and found weak or no evolution. This result can be 
interpreted in two ways. First, given the change of ~1 
magnitude of surface brightness dimming in the same 
population, the result is consistent with a passively 
evolving stellar population at a given mass, with no 
growth of galaxy disks. Second, it could imply that 
galaxies are growing, but that they must grow in such 
a way as to evolve along the stellar mass - size relation 
with time. 

On the other hand, galaxy disk formation theory pre- 
dicts that the sizes of disks, both individual and as a 
population, should be growing since z=l. In the stan- 
dard picture, gas in a halo conserves its specific angular 
momentum (equal to that of the dark matter) as it cools 
to form a centri fugally supported disk that grows from 
the i nside out ([White fc Reed 119781: iFall fc Efstathioul 
1980). In the simplest model, in which the density pro- 
file of galaxies is modeled as a singular isothermal sphere 
(SIS), the radius of the resulting disk is proportional 
to the parent halo virial radius, which grows inversely 
with the Hubble p arameter, ff(z)"" 1 . For a co ncordance 
ACDM cosmology (|Komatsu et al.ll2008ll2010h . this rela- 
tion predicts that disks at z=l should be nearly a factor 
of two smaller than their z=0 counterparts at the same 



circular velocity, V c (|Mo et al.l 119981: iMao et al.l 119981: 
Ivan den BoschHl998j). 

iSomerville et al.l (|2008f ) adopted more reasonable as- 
sumptions to update the simple theoretical model that 
suggests that galaxy disks should grow by nearly a 
factor of two since z=l, primarily by adopting NFW 
density profiles rather than SIS profiles. Unlike the 
SIS model, the concentrations of NFW profiles were 
lower in the past. Combined w ith disk stability ar- 
guments ([Efstathiou et al.l 119821 ). this model predicts 
a much weaker evolution in disk sizes ( 15%-20%) at 
a fixe d stellar mass, consistent with the IBarden et al.l 
(2005) results. These semi-analytic models (SAMs) do 
not follow the evolution of individual disk galaxies with 
time. Rather, they examine an instantaneous popu- 
lation of disks. Clearly, this is the type of informa- 
tion derived from observations, but it prevents an in- 
terpretation of the evolution in individual disk galax- 
ies based on observations of t he population as a whole. 
iFirmani fc Avila^R ccsc (2009]) instead follow the evolu- 
tion of individual disks in their SAM, and demonstrate 
that disks tend to grow along the stellar mass - size re- 
lation, with only wea k evolu tion, consistent with the re- 
sults of IBarden etHI (pOOl . 
The analytic models d iscuss e d above (e.g.. IMo et al 



19981: ISomerville et al.l 120081: IFirmani fc Avila-Reese 



20091) have until recently been the only available the- 



oretical tool with which to invest igate the evolutio n 
of galactic disk sizes (though see iBrook et al.l [2006). 
Simulations, with their ability to capture complex gas 
processes in mergers and subsequent SF, should be an 
ideal tool that allow a better trace of the distribution 
of stellar light, while SAMs must assume that light 
follows mass. However, simulations of disk galaxy 
formation in a CDM context have historically produced 
unrealistic disks that are too compact (dense) , too small 
overall, and rotating too f ast at a given radius (e.g. , 
Steinmetz fc Nayarrol 119991: iNayarro fc Steinmetel 120001: 



Eke et al.l 120011: lAbadi et all 120031: IGovernato et al l 
20041) . This failure has been named the "angular 
momentum catastrophe." This catastrophe has been 
largely associated with the overcooling problem, in 
which baryons cool too quickly at early times and 
become very dense and concentrated at the center of 
halos before merging. These halos then experience dy- 
namical friction in subsequent mergers, and the resulting 
disks show t he classic signs o f the angular momentum 
catas t rophe (|Nayarro fc Bend 119911: iNavarrq fc White! 
[1991 Katz et al.l 119941: iMaller fc Dekell l2002r K Thus, 
feedback mechanisms at early times have historically 
been in voked to prevent overcooling with moderate 
success (iDekel fc Silkl 119861: iSom mer-Lar sen et al.l 119991 : 



Thacker fc Cou chman 2000, 2001; Somm er-Larsen et al 
20031: iRobertson et al.1 120041: lokamoto et al.1 12005 



D'Onghia et al.l 120061: iKeres et all l2009t ISales et all 
20101) but additional artificial exchange of angular 
momentum can still occur between baryons and 
dark matter in SPH simulations with l ow numer- 
ical r esolution (IGovernato et al.l 120041 iNaab et all 



20071: IGovernato et alT 
Maver et al.l 120081: 



2007; 
iPiontek 



Kaufmann et al.l |2007 
Steinmetzl 120091: 
to avoid the angular 
momentum catastrophe and produce realistic disks, 
simulations must have both a physically motivated 



Sales et al.l 120101) Hence 
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Fig. 1. — Galaxies h277 (top) and h799 (bottom), seen edge on 
(left) and face on (right) at z=0 in the B band. These galaxies 
span the full range of virial masses presented in this paper, two 
orders of magnitude. 



feedback prescription, and very high numerical (mass 
and force) resolution dScannapieco et al.1 I2008L 120091 : 
iZavala et al1l200l ICeverino fc Klvpinll2009h . 

The simulations used in this study are part of a suite 
of very high numerical resolution disk galaxy simulations 
that incorporate a star formation (SF) and supernova 
(SN) feedback scheme that has been shown to overcome 
many of the past problems of disk galaxy simulations to 
successfully match a number of observed properties of 
galaxies (as discussed further below, and shown by the 
results of this paper). 

In this paper, we use simulations of individual disk 
galaxies to follow their evolution in the magnitude - size 
plane and the stellar mass - size plane. We attempt to 
mimic the observations of these galaxies in a manner sim- 
ilar to the observed population, generating mock surface 
brightness images to measure disk scale lengths, avoiding 
the assumption that light must follow mass. We demon- 
strate that the simulated galaxies have properties similar 
to real galaxies. In §2 we describe our small sample of 
very high resolution simulated disk galaxies that span a 
representative range of masses, merger histories, and spin 
values. We show that our limited sample is in excellent 
agreement with the population of galaxies surveyed at 
both z=0 (§3) and z—1 (§4). In §5, we disentangle the 
evolution of individual galaxies from the evolution seen 
in the population as a whole. We conclude in §6. 

2. SIMULATIONS & ANALYSIS 

2.1. The Simulations 

These simulations were run with the N-Body + 
Smoo thed Particle Hydrodynamics (SP H) code Gaso- 
line (jStadell [200lt IWadslev et all 12001 in a fully cos- 
mological ACDM context using WMAP year 3 parame- 



ter^ (On = 0.26, A=0.74, h = 0.73, cr 8 =0.77, n=0.96). 
The galaxies were originally selected from either a 25Mpc 
or 50Mpc (depending on their mass) N-Body simula- 
tion with uniform mass resolution throughout, and then 
resimulated at higher resolution wi th baryons using the 
volume renormalization technique (jKatz fe White! Il 9931 : 
iNavarro fe White! fl994) . This technique allows for sig- 
nificantly higher resolution on the central galaxy while 
also capturing the effect of large scale torques that are 
thought to deliver angular momentum to the galaxy 
(|Whitdll984 iBarnes fe Efstathioulll987t ). 

The SF/SN scheme used in these simulations has been 
shown to effectively regulate star formation efficiency as 
a function of halo mass, resulting in a stellar mass - 
metallicity relationship for the simulated galaxies that 
is in excellent agreement with ob s ervations both locall y 
and at high z ([Brooks et al.ll2007t iMaiolino et~aLll2008f t. 
This regulation of star formation also leads to realistic 
trends in gas fractions, with our lowest mass galaxies 
being the most gas rich (|Brooks et al.ll2007l ). reproducing 
the observed incidence ra te of Damped Lyman a systems 
(QSO-DLAs) at z = 3 ([Pontzen et all [2008). and the 
colum n densities of both QSO-D LAs and GRB-DLAs at 
z = 3 (jPontzen et al.ll2008l [20Toh . Of critical importance 
for the present study is the fact that these simulated 
disks maintain sufficient angular mom entum to match 
the observe d Tully-Fisher relationship ([Governato et al.l 
2008, 2009) and produce galaxies with realistic disk sizes, 
as shown below. 

The full details of our physically motivated SN 
feedback imp l emen tation were originally presented in 
iStinson et al.l ([2006f ). Briefly, the SF prescription en- 
sures that the SFR density is a function of gas den- 
sity according to the observed slope of the Kennicutt- 
Schmidt law, and a SF efficiency parameter, c*, sets the 
normalization of this relation. Each star particle repre- 
sents a simple stellar population, born w ith a Kroupa 
initial mass function ([Kroupa et al.l fl993h . As massive 
stars go SN, energy and metals are deposited into the 
nearest neighbor gas particles. The SN feedback recipe 
calculates the radius affected, and turns off cooling in 
those affected neighboring gas particles until the end of 
the sno wplow phase as described by the Sedov- Taylor so- 
lution ([McKee fc OstrikeHll977T) . The amount of energy 
deposited amongst those neighbors is 0.4xl0 51 ergs, as 
was adopted in all of our previous work mentioned above. 
Additionally, we include a uniform UV cosmic back- 
groun d following an updated model of lHaardt fc Madaul 

ii. 

iGovernato et ahl (|2010f ) demonstrated that when force 
resolutions <100pc can be achieved, high density 
peaks that mimic SF regions in giant molecular clouds 
can be resolved with several hund reds of gas parti- 
cles in ~10 6 M(n clump s (see also iBooth et all (20071: 
ICeverino fc Klypinl 120091 ) . This allows for the adop- 
tion of a realistic density threshold for SF (100 
amu/cm , iRobertson fc Kravtsovll2008t iTasker fc Brvanl 
120081 iSaitoh et al.l I2008D . To match the Kennicutt- 
Schmidt law, the high density SF threshold must be off- 
set by a slightly higher value of c* = 0.1. This pre- 
scription leads to enhanced gas outflows that remove 

11 The choice of WMAP3 cosmology over WMAP7 cosmology 
has no impact on the results presented here. 
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low angular m omentum gas from th e central regions 
of the galaxy (|Governato et al.l 120101 Brook et al., in 
prep.). At low masses, this creates a bulgeless disk with 
a linearly rising rotation cur ve, comparable to those ob- 
served (Oh et al., in prep; Ivan den Bosch et~aH [20011 
Ide Blok et al~l l2008). The low mass galaxies used in this 
paper, drawn from a 25Mpc volume, can achieve force 
resolutions ^lOOpc. Hence, these simulations adopt the 
more realistic feedback prescription with high density 
SF and c* = 0.1. However, the more massive galaxies 
presented here are drawn from 50Mpc volumes, making 
it computational expensive to achieve similar force res- 
olutions in a reasonable time. Hence, these simulated 
galaxies cannot resolve the high density SF peaks easily, 
and a lower SF density threshold must be adopted, 0.1 
amu/cm 3 , with a lower c* = 0.05. These latter values 
have been ad opted in all of our previous work on MW 
mass galaxies (iBrqoks et al.ll2007l: iGovernato et al.ll2007l: 



Brooks et al. 2009; G overnato et al. 2009k iPontzen^ et al.l 



20081 12010ft . and are a compromise that allow realistic 



disks to form, but are inefficient at driving gas outflows. 
Thus, low angular momentum material that might be 
lost from the central regions is maintained, helping lead 
to the creatio n of large bulges in the MW mass galax- 



20011 Ivan den Bosch et al 


2002; Mailer & Dekell 


2002; 


van den Bosch et al.l [2003 
Duttonl 1200*917" Additional 


iD'Onghia & Burkertl 
y, the creation of these 


2004; 
large 



bulges can be due to missing physics. In particular, AGN 
feedback has not yet been added to these simulations, 
and is potentially a key mechanism to help create real- 
istic bulges. Due to our small sample size, it is difficult 
to quantify if our bulges are substantially different from 
observed bulges in the same galaxy mass range. 

The simulated galaxies used in this work were selected 
to span a range of spin values and merger histories, with 
the last major merger redshift as low at 0.8. Properties 
of these simulations are listed in Table [TJ 

2.2. Deriving Disk Scale Lengths 

All of our simulations are of field galaxies that do not 
undergo major mergers below z ~0.5, making them disk 
dominated at redshift zero. In order to compare our 
simulated disk scale lengths to observational results in 
as realistic a way as possible, we wish to fit the light 
profile rather than the underlying mass profile. We cre- 
ate artificial surface brigh tness images o f our simulated 
galax ies using Sunrise (jJonssonl 120061 : Uonsson et al.l 
2010), a Monte Carlo radiative transfer program that 
produces a spectral energy distribution (SED) for each 
resolution element of the simulations. This is done by 
identifying the age and metallicity of each star parti- 
cle, which are then convolved wit h the Starburst99 stel- 
lar population synthesis m odels (|Leitherer et al.l 119991 : 
iVazquez fc Leithererll2005|) to produce an SED. SUNRISE 
assumes dust tracks with the metallicity of the gas par- 
ticles, and performs ray tracing from each star particle 
to compute the observed SED including absorption and 
scattering. From this, we generate mock images in cho- 
sen filter bands. Figure [TJ shows artificial edge-on and 
face-on B band images, including dust reprocessing, for 
two different mass galaxies. 

Most of the observational surveys that have investi- 
gated the evolution of disk sizes worked in the rest frame 



B band across redshifts. To compare to these results as 
a function of redshift, we have generated face on i band 
surface brightness images at z=0.5 and z=l that include 
the effects of surface brightness dimming, and using rest 
frame face on B band surface brightness images for the 
z=0 galaxies. Two of our disk galaxies are undergoing a 
major, disruptive merger at z—1. For these galaxies, we 
created face on images for the most massive progenitor 
at a time just prior to the merger (z=1.25) while the disk 
was still rather undisturbed. 

The face on Sunrise images were fit with a Sersic 
bulge component and an exponential disk component us- 
ing the publi cly available 2D surface brightness fitting 
code Galfit ()Peng et al.ll2002f ). Additionally, bars exist 
in both h603 and h986 at z=l that were fit by a third 
component, reducing the B/D ratio of these galaxies at 
this step. The Galfit results were checked against a ID, 
radially averaged fit generated from the same images, and 
found to be in good agreement. To mimic observations, 
the fits were required to be a good match down to a lim- 
iting surface brightness value of ~26 mag/arcsec 2 (hence 
the importance of using surface brightness dimmed im- 
ages at z=l). Example fits are shown in Fig. [2] for the 
same galaxies shown in Fig. Q] with K band fits shown 
simply for a comparison of the B band results to those at 
a longer wavelength. We verified that using redshifted, 
surface brightness dimmed Sunrise images in the i band 
at higher z yielded similar disk scale lengths as the rest 
frame B band. 

As can be seen in Table [IJ a couple of these galax- 
ies have undergone major mergers since z=l, while some 
galaxies have had a very quiescent history. We note that 
at z—1, two of our galaxies (h239 and h285) had no com- 
ponent that could be well fit by an exponential. This 
may be due to mergers that occur in these galaxies near 
this time. Galaxy h239 has the most active merger his- 
tory of our simulated galaxies, with continual mergers 
(both major and minor) until z ^0.5. Galaxy h285 be- 
gins to accrete a satellite with 1/7 of the stellar mass 
of the main halo at z=1.25. However, the core of this 
satellite does not merge with the main halo until z ~0.8. 
In both cases, we searched in 250 Myr periods from 0.75 
< z < 1.25 to identify a step with a possible exponential 
disk, but none could be found. Hence, those steps have 
b een excluded from t he hig h redshift analysis below. 

iScannapieco et aTl (|2010f ) demonstrated that measur- 
ing disk scale lengths based on light rather than a kine- 
matic decomposition could dramatically increase the re- 
sulting disk-to-total, D/T, ratio of their simulated galax- 
ies, making measurement technique a potentially impor- 
tant reconciler between observations and simulations. In 
Fig. [3] we plot the z=0 surface mass density as a func- 
tion of radius for both the stars (solid line) and the gas 
(dashed line) in the simulated galaxies. Also shown as 
the dotted line is the exponential fit derived in the i band 
for these simulations (listed in Table [T]). In a future pa- 
per (McCleary et al., in prep), we will explore in detail 
the differences between mass/kinematic results and light 
results (e.g., scale lengths, bulge-to-disk ratios, the role 
of dust and inclination, etc.). Here, however, we note 
that an initial comparison of mass versus light exponen- 
tial scale lengths demonst rates a less dramati c differ ence 
in our simulations than in IScannapieco et aTl (|2010D . As 
seen in Fig. [3l the i band scale length is generally a good 
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TABLE 1 

Simulated Galaxy Properties 



simulation 


M vir 


M» 


M DM 
particle 


M sph t - , 

parttcle 


A 


A 9 


z imm 


e 


N within K v i r 




M 


M 


M© 


M 








pc- 


dm+star+gas 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


CS) 


(10) 


h516 


3.9xl0 10 


2.6xl0 8 


1.6xl0 4 


3.3xl0 3 


0.05 


0.05 


1.2 


87 


3.5xl0 6 


h799 


2.2xl0 10 


1.3x10 s 


1.6xl0 4 


3.3xl0 3 


0.04 


0.05 


3.0 


87 


1.9xl0 6 


h603 


3.8X10 11 


3.1xl0 10 


3.0xl0 5 


6.3xl0 4 


0.06 


0.11 


1.1 


231 


3.8x10 s 


h986 


2.1X10 11 


2.2xl0 10 


3.0xl0 5 


6.3xl0 4 


0.04 


0.07 


0.8 


231 


2.4x10 s 


h239 


8.5X10 11 


7.8X10 10 


1.2xl0 6 


2.1x10 s 


0.03 


0.05 


1.1 


347 


2.8x10 s 


h258 


8.0X10 11 


7.4xl0 10 


1.2xl0 6 


2.1x10 s 


0.04 


0.07 


0.8 


347 


2.8x10 s 


h277 


7.1X10 11 


6.9X10 10 


1.2xl0 6 


2.1x10 s 


0.03 


0.04 


3.0 


347 


2.3x10 s 


h285 


8.7X10 11 


8.1xl0 10 


1.2xl0 6 


2.1x10 s 


0.02 


0.05 


1.9 


347 


3.0x10 s 



Note. — Properties of the galaxies as drawn from the simulations. Columns (2) and (3) list the virial 
mass and total stellar mass of the halos at z=0. Columns (4) and (5) list the mass resolution of individual 
dark matter and st ar particles, respectively. Column (6), A, is the dimensionless spin parameter, a la 
Bullock ct al. (2001), for the entire halo. Column (7) lists the spin parameter for all gas within the halo. 
The last major merger (1mm) redshift in column (8) is defined at the time when the cores merge of two 
galaxies initially ~3:1 in halo mass. Column (9), e, is the spline gravitational force softening. The final 
column (10) lists the total number of particles within the virial radius of the halo at 2=0. 
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Fig. 2. — Face-on, radially averaged B band (left panels) and K band (right panels) surface brightness profiles (solid lines) at z=0 for 
the same galaxies as in Fig. [T](h277 and h799). Bulge/disk decompositions resulting from Galfit are shown, though the low mass galaxy 
in the bottom panel is bulgeless. The long dashed line in each panel at 26 Ms/arcsec 2 represents the limiting surface brightness out to 
which the profiles were fit. 
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Fig. 3. — The z=0 surface mass density distribution for the sim- 
ulated galaxies. Stellar surface density is shown by the solid lines, 
while gas surface density is shown by the dashed lines. For com- 
parison, the dotted lines in each panel show the i band scale length 
derived using Sunrise/Galfit for each galaxy (listed in Table|5]|. 

match to the stellar surface mass density. For a quantita- 
tive comparison, scale lengths were derived for the stellar 
mass over the same region that the disk dominates the 
light profile. The two methods were found to agree to 
within 10% for 75% of the galaxies. Variation in the 
techniques appears to be mostly attributable to varying 
mass-to-light ratios (McCleary et al., in prep). 

The bulge/disk decompositions were performed on the 
dust extinguished Sunrise results. However, disk scale 
lengths and central surface brightnesses were found to be 
identical in both the face on Sunrise images with and 
without dust reprocessing, while the central bulge com- 
ponent was dust extinguished. Thus, using the dust ex- 
tinguished versus non dust extinguished decomposition 
has no effect on the results for the disk analysis pre- 
sented below. It is generally expected that dust should 
have some effect on both the central surface brightness 
and the scale length of the disk c o mponent, even for 
face on orientations (jCalzettil 120011 : iGraham fe; Worleyl 
2008). We speculate that our more massive galaxies 
have little/no effect from dust due to the fact that they 
are also more gas poor than observed galaxies (because 
they use up gas forming too many stars, discussed in 
more detail in Section I5T2]) . Our low mass galaxies are 
gas rich, in agreement with observed galaxies in thi s 
mass range dGovernato et al.l 



20101: IGeha et al.l [2006), 



2007D . making dust neg- 



but metal poor (|Brooks et al.l 
ligible. 

It is important to note that the following analysis 



purposely focuses on the evolution of the disk alone. 
That is, the discussion below of how properties evolve 
with time neglects the bulge component of the galaxy. 
This is done for three reasons. First, as discussed 
above, our MW-mass galaxies tend to have larger bulges 
than observed. This is a problem that ha s historically 
plagued MW-mass simulated disks (e.g., lAbadi et al.l 
l2003t IScannapieco et aLll2009D . and is likely due to the 
inability to correctly model SF and feedback at the nec- 
essary high resolutions to drive loss of low angular mo- 
mentum gas, and due to missing physics such as AGN 
feedback. Hence, while the disks of our simulated galax- 
ies appear to be in good agreement with observed disks 
(as presented below), the growth of the bulge in these 
massive galaxies is probably not modeled correctly. Sec- 
ond, observations of the evolution of disk galaxies in the 
magnitude - size plane have either a) selected galaxies 
that can be fit with a small pure Sersic index, n, so that 
they are disk dominated, or b) done a bulge/disk decom- 
position. In either case, the results we wish to compare 
to are concerned with disks only. Finally, our goal here is 
to examine the growth of disks using these simulations, 
not bulges or spheroidal components. 

"Observable" properties of the simulated galaxies are 
listed in Table [2J Total magnitudes in each band are de- 
rived from Sunrise, with each galaxy being integrated 
out to roughly 50 comoving kpc (i.e., essentially all of 
the flux from the simulated galaxy is contained). Ta- 
ble also lists the z—0 stellar mass derived for the en- 
tire galaxy (disk + spheroid) based on Sunrise colors 
usin g the k_sdss_bell routin e in the KCORRECT pack- 
age (|Blanton fc Row cis 2007). A comparison to Table Q] 
shows that the simulated stellar masses are generally 
^40% larger than those derived from the photometric 
results. A more detailed comparison will be presented 
in McCleary et al., (in prep). In Table [2] we list the 
bulge-to-disk ratios (B/D) despite the caveat mentioned 
above regarding large bulges in our higher mass galax- 
ies. Disk scale lengths in the i band at z=0 are given 
in Table [21 while the B band scale lengths used in the 
following plots are given in T able |3l As seen in observa - 
tional studies (|de Jongl Il996al iMacArthur et al.ll2003f l. 
the B band scale lengths tend to be longer than the 
i band scale lengths, at least for the six most massive 
galaxies (with the two dwarf galaxies exhibiting different 
behavior). It has been debated whether this trend for 
disk scale lengths to be longer at sh orter wavelengt hs is 
due to dust effects or age gradients (jde Jong||1996bD . As 
noted above, dust plays no role in our face on disk fits, 
meaning that this trend in the simulations is due entirely 
to age and metallicity gradients. 

3. THE FUNDAMENTAL PLANE FOR DISK GALAXIES 

Size, magnitude (or mass), and rotational veloc- 
ity make up a fundamental plane for disk galaxies 
dPizagno et alj|2005t iGnedin et al.ll2007HCourteau et al.l 
2007t). iGovernato et all (|2009l see their figure 5) demon- 
strated that these simulated galaxies lie on the observed 
magnitude - velocity relation for disk galaxies (also 
known as the Tully-Fisher relation). In this paper we 
derive disk scale lengths in order to extend this analysis, 
and show that our galaxies are also a good match to the 
observed magnitude - size and velocity - size relations 
as well. All of our galaxies have sizes and velocities (Ta- 
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TABLE 2 

"Observable" Galaxy Properties 



simulation 




M s 


B/D, 


B/D B 


g-r 


ML orr 




V 2 .2 














M 


kpc 


km/s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


h516 


-16.8 


-15.9 


0.08 


0.14 


0.53 


3.4xl0 8 


1.0 


51 


h799 


-16.2 


-15.6 


0.00 


0.00 


0.37 


1.1x10 s 


0.7 


43 


h603 


-21.3 


-20.5 


1.25 


0.72 


0.52 


2.2xl0 10 


1.0 


113 


h986 


-21.1 


-20.3 


0.63 


0.38 


0.45 


1.3xl0 10 


2.8 


137 


h239 


-22.5 


-21.8 


0.35 


0.16 


0.41 


4.2xl0 10 


2.5 


246 


h258 


-22.4 


-21.5 


0.86 


0.45 


0.50 


5.3xl0 10 


5.3 


204 


h277 


-22.1 


-21.2 


0.63 


0.42 


0.52 


4.3xl0 10 


2.3 


250 


h285 


-22.3 


-21.3 


1.33 


0.78 


0.55 


5.6X10 10 


4.5 


203 



Note. — Magnitude and color results are dust free measurements for 
the entire galaxy (disk and spheroidal components) at 2=0. B/D ratios in- 
clude the effect of dust extinction (see discussion in Section |2 . 2D . The stellar 
mass es listed in column (7) a re derived using k_sdss_bell in the kcorrect pack- 
age (Blanton & Rowcis 2007) from SDSS ugriz colors provided by Sunrise. 
Rotation curve velocities are measured at 2.2 i band disk scale lengths. 



bleEJ that agree well with the observed velocity - size 
trend (|Courteau et al.l [20071. Two of the more massive 
galaxies have slightly high velocities for their size (h239 
and h277), but still lie within the 2<r scatter of observed 
galaxies. 

The B band disk magnitudes, Mb, and scale lengths, 
hg, of our simulated galaxies at z=0 are shown in 
Fig. |U The simulated galaxies (red circles) are com- 
ared to three observational sa mples, iMacArthur et al.l 
2003) an d IGraham fc Worlevl ((20081) at the massive 
end, and Ivan Zed (120001) at the low mass end. The 
IGraham fc Worlevl ()2008| ) data have adopted a dust 



correction based on [ Driver et al.l (|2008[ ). while the 
IMacArthur et al.1 ()2003l ) data are dust corrected based 
on inclination, foll o wing equations 1, 2, and 7 listed in 
IGraham fc Worlevl (|2008f ). Again, possibly due to low 
gas (and therefore dust) content, the six more massive 
galaxies shown here are free of dust effects. As such, they 
should be compared to dust cor rected observat ional data. 
At the low mass end, however. Ivan Zeel (f20M ) makes no 
internal dust corrections, as the low metallicities of the 
dwarf galaxies in her sample are expected to lead to little 
effect from dust Q Clearly, the disk sizes of the simulated 

12 Note that there appears to be a slight offset between the dwarf 
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galaxies are in good agreement with observed disk sizes. 

We use the B band scale lengths of the observational 
samples in Fig. [4j which can be compared to the i band 
results at z=l in the next section since the B band is 
rcdshifted into the i band at z=l. However, the mag- 
nitude - size (and stellar mass - size) relation has been 
measure d in multiple bands for a large sam ple of SDSS 
galaxies (IShen et al J 1 20031: iFathi et alJl2010h . The SDSS 
results are consistent with the results plotted in Fig. [4] 
though we purposely select observational samples that 
have done a bulge/disk decomposition and fit a pure ex- 
ponential to the disk surface brightness fit, in order to 
isolate the evolution of the disk and eliminate any con- 
tamination from a central spheroid. 

As discussed in the Introduction, historically it has 
been a challenge for cosmological disk galaxy simulations 
to match t he observed fundamental rel ations for galaxy 
disks (e.g.. iNavarro fc Steinmetl l2000f ) . The failure to 
reproduce observed trends is the result of the angular 
momentum catastrophe, exacerbated by the fact that 
previous studies followed mass/kinematic results rather 
than light profiles. 

Our success at matching the observed disk scaling re- 
lations is a result of the fact that our disk baryons 
maintain angular momentum through time, due to an 
increased resolution a nd a physicall y motivated feed- 
back mechanism (see Brooks ~ 2010, for a review of 
how each of these processes contribute to the forma- 
tion of realistic disks). A number of works have at- 
tempted to iso late the separate role of resolution and 
feedback (e.g.. | Governato et al.l l2004t IKaufmann et al.l 
120071: iNaab et all 120071; iPiontek fe Steinmetzl l2009f ) . "hi 
N-Body + SPH simulations, dark matter particles are 
typically an order of magnitude more massive than the 
gas particles, leading to an exchange of kinetic energy in 
two body interactions that kinematically heats the disk, 
randomizing velocity vectors, and potentially turning 
a disk component into a spheroid ( Stci nmetz fc Whitel 
119971; IMaver et all 120081) . This effect can be reduced 
at higher mass r esolutions that lower th e discrepancy in 
particle masses. IKaufmann et al.l (|2007f ) used controlled 
(non-cosmological) N-Body+SPH simulations to evalu- 
ate the amount of angular momentum that can be lost 
purely due to resolution. They concluded that 10 6 par- 
ticles within the virial radius are necessary for a disk 
galaxy to maintain roughly 90% of its original angular 
momentum. 

Yet even if resolution can be increased to the point of 
maximum angular momentum retention, disks will still 
suffer dramatic angular momentum loss compared to ob- 
served disk galaxies if energy feedback is neglected. Feed- 
back prevents rapid, early cooling of gas particles (the 
overcooling problem discussed in the Introduction; see 
references therein). Heating and expansion of the gas 
creates a hot reservoir, allowing the gas to cool at later 
times after the era of rapid mergers, and preventing an- 
gular inomfintum_l_oss |_via_ dynamical friction in mergers 
(e.g.. Mailer fc Dekelll2002l) . Additionally, by preventing 
overcooling, feedback prevents gas from rapidly turning 
into collisionless star particles. Without feedback, ha- 



galaxy sample and the higher mass sample. This break has been 
observed previously, and sugge sts a structural d ifference between 
dwarf and larger disk galaxies ( Schombcrt 2006). 



los of all masses are equally efficient at converting gas 
into stars (jBrooks et all 120071) . producing galaxies that 
are too gas poor compared to z — disk galaxies. Be- 
cause feedback regulates star formation, it allows for gas 
reservoirs to develop that allow disks to survive to the 
present da y (Hopkiis et al. 20091: iGovernato et al.l 120091 : 
iMoster et al.ll2010aD . The existence of large, thin disks 
at the present day thus requires feedback. 

Not all SN feedback schemes lead to disk galaxies 
that satisfy observed constraints such as the Tully- 
Fisher relation, size - velocity or size - luminosity 
relations, or the stellar mass - metallicity relation. 
Because the nearest neighbor gas particles surround- 
ing SNe are dense and cold, a simple energy depo- 
sition will quickly b e radiated away and not affect 
the s imulation (e.g., iKatd 1992; S teinmetz fc Navarro! 
fl999h . Two main schemes have been adopted to 
overcome this problem. In the first, a multiphase 
model of the ISM is impleme nted with a hot gas 
reservoir inside each gas particle dHultman fc Pharasvnl 
1999t iMarri fc White! I2O03I; iSpringel fc Hernquistl 12005 



Scannapieco et al.l 120061 : lHarfst et alll2006l ). overcoming 



the problem of thermal energy being spread over the 
entire particle mass, and preventing the hot gas parti- 
cles from being artificially influenced by their cold gas 
nearest neighbors. In the second, cooling is turned off 
in the gas particles near a SN explosion in order to 
mimic the sub- resolution adi a batic expansion of the SN 
(jThacker fc Couchmanl 120001 . 120011 ). The "blastwave" 
scheme adopted here also turns off cooling in nearest 
neighbor particles, but attempts to model this based on 
what is known about actual SNe, determining the ra- 
dius of each SN remnant base d on the analytic blast- 
wave solution for a S N remnant (iMcKee fc Ostrik"er1ll977l : 
iStinson et all 12006). and cooling is only turned off for 
those particles within the blast radius. Because many 
SNe typically contribute feedback within a dense star 
forming region, the thermal energy from all of these 
SNe can combine to create a larger blast radius. The 
differences in the resulting disk between a simple en- 
ergy deposition ( "thermal" fee dback) and the blast wave 
model have been exam ined in IMaver et al.l (|2008l ) and 
IGovernato et al.l (|2008). As examined by these previous 
works, the adoption of the "blastwave" feedback model, 
combined with high resolution (all simulations presented 
here have more than 10 6 particles within the virial ra- 
dius at z — 0, see Table [T]), overcome past problems with 
cosmological disk galaxy simulations and allow for the 
present study. 

4. EVOLUTION OF THE MAGNITUDE - SIZE RELATION 

Having established that these simulated disk galaxies 
match observed disk scaling relations at z—0, we now 
use them to investigate the degeneracy in the magnitude 
- size evolution with time. First, we establish that the 
simulations also match the available observational data 
at z=l. 

4.1. Evolution as a Population 

The evolution of the magnitude - size relation for these 
simulated galaxies is shown in Fig. [5J The top panel 
shows the simulation results at z=0, 0.5, and 1. It is evi- 
dent that the magnitude - size relation for these galaxies 
is evolving in time, dimming in surface brightness since 
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Fig. 5. — Rest frame B band disk scale length as a function of 
magnitude for our simulated galaxies. Top panel: The evolution of 
the population with time. Middle panel: Simulated points at 2=1 
as in the top panel, but now compared to observ ational results 
for decomposed disks from MacArthur ct al. (2008) and Miller et 
al. {in prep). Bottom panel: Same as Fig. [4] shown again for com- 
parison to the high z results. 



z=l. Despite the fact that these are individual galaxies 
being followed in time, the snapshot of the population 
of disks at each redshift is consistent with the observa- 
tions. This is evident in the bottom two panels. The 
middle panel compares the simulated disk scale lengths 
at z—1 to the decomposed i b and disk scale len g ths o f 
galaxies at 0.9 < z < 1.2 from iMacArthur et all (|2008f ) 
and Miller et al. (in prep). No dust corrections are ap- 
plied to the observational z=l data, as these galaxies 
are expected to be low metal licity where correctio ns are 
uncertain and often ignored (jDutton et a l. 2010b). Any 
realistic dust corrections are likely to be tiny, and not 
affect the results. The bottom panel is a reproduction 
of Fig. |4l shown again for easy comparison of the results 
across redshifts. While observational biases allow for a 
direct comparison only at the massive end at z—1, the 
disk sizes of the simulated galaxies are in good agree- 
ment with observed disk sizes, as a function of redshift. 
We note that observatio nal data that used a single Ser- 
sic component fit (e.g., iBarden et al.l 120051 ) are also in 
good agreement with the z=l data presented here (and 
lie in the same magnitude range), despite the lack of a 
bulge/disk decomposition. For clarity, we show only the 
decomposed data in Fig. [SJ 

The top panel of Fig. [5] shows that there has been ap- 
parent dimming of the surface brightness of the simulated 
population since z=l. Fig. [6] quantifies this evolution, 




Fig. 6. — The evolution of surface brightness as a function of 
redshift. The two lowest mass galaxies are shown separately from 
the high mass galaxies. Data points at each z show the mean of 
the population, and the error bars reflect the standard deviation 
of the sample. 

showing the change in the surface brightness, /is, of the 
simulated galaxy disks. We have defined [Ib such that 



1-iB = M B 



disk 



5logh B + 2.Uog(2n). 



(1) 



Data points in Fig. [5] at z=0, 0.5 and 1 are for the 
mean [Ib for the population at that time, and error 
bars reflect the standard deviation of the sample. Be- 
cause our two dwarf galaxies are too faint to be gen- 
erally observable at z=l, Fig. [5] shows the evolution 
of these two galaxies separately from the higher mass 
galaxies. However, the evolution is roughly similar in 
both mass ranges, with about 1.5 magnitudes in sur- 
face brightness dimming between z—1 and z—0. A 
dimming of 1.5 magnitudes is consistent with observa- 
tional studies that have examined the surf ace brightness 
evolution of disk dominated galaxies (e.g.. iSimard et all 
1999HRavindranath et alll2004l iTruiillo fe Aguerrill2004l: 
Barden et al.ll2005t Melbourne et alj|2007t iKanwar et all 
20081) . 



4.2. Evolution of Individual Galaxies 

The previous section examined the population of sim- 
ulated disks in terms of surface brightness evolution. We 
now decompose the surface brightness trends into lumi- 
nosity and size evolution, and demonstrate how individ- 
ual galaxies evolve in each property with time. Table [3] 
quantifies the changes in the disk B band central surface 
brightness (hb,o, as opposed to (ib plotted in Fig. |6j), 
size, and magnitude for each of our simulated galaxies. 

4.2.1. Luminosity Dimming since z = 1 

It is widely expected that galaxies should undergo dim- 
ming since z—1, due to declining star formation rates. 
The key question is how much galaxies dim. It is al- 
ready apparent from Fig. [5] that our more massive galax- 
ies are undergoing growth in size between z=l and z=0, 
meaning that the change in surface brightness of ~1.5 
magnitudes over this time cannot be entirely due to lu- 
minosity dimming at a fixed size. The changes in size 
and magnitude that are quantified in Table [3] are also 
visualized in Fig. [7] Now it can be seen that there is in 
fact a dependency on mass in the evolution of the galax- 
ies. Our more massive galaxies are growing in size and 
generally undergoing only a small amount of dimming 
in total magnitude; some are even getting brighter with 
time. Meanwhile, there is little change in the scale length 
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TABLE 3 

Evolution of Galaxy Disk Properties 



simulation 


A*B 




M s 


AM B 




Ah s 


M* 


AM, 




mag/LT 


Amag/D" 


mag 


Amag 


kpc 


2=0/2=1 


M 


2 = 0/2 = 1 




2=0 


to 2 = 1 


2 = 


to 2=1 


2=0 




2 = 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


h516 


22.7 


1.5 


-15.7 


2.0 


0.95 


1.0 


4.3x10 s 


1.69 


h799 


22.1 


0.7 


-15.6 


0.0 


0.6 


1.3 


1.8x10 s 


1.32 


h603 


22.5 


2.8 


-19.8 


1.0 


1.2 


2.4 


1.3xl0 10 


2.78 


h986 


21.2 


1.0 


-19.9 


0.2 


3.0 


1.9 


l.OxlO 10 


3.22 


h239 a 


19.3 


N/A 


-21.5 


N/A 


3.2 


N/A 


9.0xl0 9 


10.17 


h258 a 


21.3 


1.6 


-21.1 


-0.2 


6.3 


2.7 


3.2xl0 10 


2.53 


h277 


19.7 


0.6 


-20.8 


0.4 


2.4 


1.4 


1.9xl0 10 


1.95 


h285 


21.4 


N/A 


-20.6 


N/A 


4.6 


N/A 


2.0xl0 10 


3.07 



Note. — For columns (3) and (5) (Afis,o an< i AMs), positive values represent dimming 
between 2=1 and 2=0. A negative value means that the galaxy disk is brighter at 2=0. Note 
that values listed in this table are for the disk only, while those in Table [2] are for the entire 
galaxy (disk and spheroid). h239 and h285 do not have exponential disks at 2=1. The stellar 
disk masses have been measured directly from the star particles in the simulations, based on a 
kinematic disk decomposition. See Section [5] for details. 
a Measured at z=1.25; undergoing major merger at z=l. 
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Fig. 7. — Left panel: Evolution in the magnitude — size plane with time. Each color connects an individual galaxy at 2=1, 0.5, and 0. 
The higher z steps are marked by diamonds, and the 2=0 step is indicated by a cross within the diamond. Note that low mass galaxies 
evolve more in luminosity (dimming), while higher mass galaxies evolve more in size. Right panel: Evolution in the stellar mass - size plane 
as a function of time. Note that these galaxies grow so that they remain on approximately the same relation with redshift. 



of our lowest mass galaxies, but their evolution in lumi- 
nosity is erratic (for this sample of two galaxies) due to 
their bursting star formation histories. 

The luminosity evolution of the low mass dwarf galax- 
ies is tied up with their particular star formation histo- 
ries. While h799 is quiescent between z=l and z=0.5 
so that it dims, it undergoes another burst of SF at low 
z, increasing its magnitude. On the other hand, h516 
stays relatively quiescent between z=l and z=0, so that 
it undergoes significant luminosity dimming (of 2.0 mag- 
nitudes). This large amount of dimming may at first 
seem to be at odds with the fact that low mass galax- 
ies are generally undergoing more SF today than they 



were in the past (this is discussed further in Section [5J 
and is one form of galaxy "downsizing"). However, this 
extreme dimming for low mass galaxies has been seen 
in observational res u lts as well ((Melbourne et al.l 120071 : 
IKanwar et all 120081 ). Melbourne et al.l (|2007t) concluded 
that small galaxies (those with half light radi£3 2 kpc) 
have undergone roughly 2.5 magnitudes of dimming since 
z=l. Given the bursty SF nature of the two low mass 
dwarfs we have in our sample, it is not possible to com- 
pare statistically to the observations, but there is noth- 

13 The half light radius is a factor 1.68 times larger than the 
disk scale length for an exponential disk. 
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Fig. 8. — The growth of the simulated disks, normalized to their 
high z sizes. The dashed and dotted lines show simple predictions 
for the growth of disk galaxy populations since 2=1. The colors 
here are the same for individual galaxies as in Fig. [7] 

ing to su ggest that our simulated dwarfs are inconsistent 
with the [Melbourne et al.l (120071) results. 

At the higher mass end. iMelbourne et al.l (|2007l ) find 
~1.5 mag of dimming since z=l for galaxies with half 
light radii > 3 kpc (disk scale lengths > 1.8 kpc), if the 
strictest interpretation is made that all surface brightness 
evolution is in luminosity rather than size. Clearly, this 
is not true for our high mass galaxies, which are growing 
with time. 

If galaxies are growing with time, then comparing 
galaxies that are a given size (e.g., 3 kpc) at z=0 with 
galaxies that are of similar size at z—1 means that similar 
galaxies are not being compared. Those z=0 galaxies will 
have been smaller at z=l, and it is unknown from our 
simulation results what size galaxies with scale lengths 
larger than ^2.5 kpc at z=l will have evolved to by z=0. 
In summary, if galaxies are growing since z=l, this will 
mimic a larger change in luminosity on the magnitude 
- size relation than actually occurs (ITruiillo fc Aguerril 
l200i IBarden et al.ll2005t iTruiillo fc Pohlenll2005D . 

4.2.2. Change in Size since z = 1 

Our simulated disk galaxies as a population appear 
to match the observational results for the evolution in 
the size-magnitude plane very well (Fig. [5]) . Yet Fig. [7] 
demonstrates what the observations cannot witness: in- 
dividual galaxies are growing with time, so that the evo- 
lution in the magnitude - size plane cannot be due to 
dimming alone. Disk scale length evolution for each 
galaxy is plotted in Fig. [8j normalized to their size at 
z=l. These values are also quantified in Table[3J Only a 
few of these galaxies show negligible size evolution, with 
the two lowest mass galaxies undergoing the smallest size 
changes. The remaining galaxies show clear growth, and 
about half evolve by the amount predicted by the simple 
SIS model, (nearly a factor of two back to z—1), or more. 
The dotted line in Fig El shows the gr owth predicted for 
the simple SIS model ()Mo et al.H"l998l ) while the dashed 
line shows the growth determined by iMao et al.l (1998) 
using 16 galaxies at 2=1. 

However, it is not necessarily expected that any given 
galaxy should follow th e Hjz) -1 analyt i c grow th pre- 
dicted by the SIS model. iSomerville et al.l (|2008| ) demon- 
strated that more sophisticated models adopting NFW 
profiles could bring the observations and theory closer 
in line. The prediction that sizes scale with H(z)^ 1 is 



based on the assumption that 

r 2 oo = V c /WH(z) (2) 

where V c is the circular velocity of the halo and T2oo is the 
radius at which the mean density is equal to 200 times 
the critical density. Then 

h cx Ar 20 o (3) 

where A is the dimensionless spin parameter ()Peeblesl 
I1969H . If the disk size follows H(z)~ l , both A and V c 
must stay roughly constant with time. While V c does 
stay roughly constant back to z—1 in these galaxies, A 
does not. It is known that A may vary fo r individual ha- 
los due to mergers or smooth accretion (iGardnerl [200ll: 
[ Mailer et all 120021: IVitvitska et al l l2002t iPeirani et al. 
2004| iHetznecker fc BurkertJ 120061: iD'Onghia & Navarre 
l2007f h The previous work of analytic models do not re- 
quire A to stay constant, as they simply use a snapshot 
in time and the spin values of halos at that snapshot. 
That is, they study a population of galaxies at a given 
time, and make no attempt to follow an individual galaxy 
with time as we do here. Hence, it is not surprising that 
some of the simulated galaxies presented here deviate 
from growing precisely as H(z) . 

A better question, if we want to compare our disk sizes 
to analytic models, is to ask if the scale lengths measured 
based on our light profiles are similar to those predicted 
by analytic models. For an NFW density distribution, 
the disk scale length, is 

R d = —}=—r W of R (\,c,m d ,jd) (4) 

where c = r2oo / r s measures the halo concentration, and 
r s is NFW the scale radius. f c uses c to determine the 
energy of the NFW halo compared to the SIS model, 
and fR(X,c,m,i,jd) is a factor that accounts for the 
gravitational effect of the disk (adiabatic contraction). 
3d = Jdisk/ J200 and m d = M d i S k/M 2 Qo are the fraction 
of disk angular momentum and mass, respectively, to the 
total halo. For sizes to scale as H(z)^ 1 , it is assumed 
that jd/md = 1. This is likely to no t be true, though, as 
demonstrated bv lSales et al.l (|2009P i. 

We have determined the values of f c , uid, and jd 
for each of our galaxy halos at z=0, directly from the 
simulation results. We compared the Rd expected from 
these quantities, as derived by equation [U to the i band 
scale lengths listed in Tabled] The results agree to within 
10% for only the two lowest mass galaxies. For the more 
massive galaxies, the predicted NFW Rd is consistently 
shorter than that found for the i band, and for four of 
the six massive galaxies, the measured i band results are 
a factor of 2-4 larger than the NFW prediction. 

The main source of this discrepancy appears to be 
in the value obtained for the spin value, A, for the en- 
tire halo. If the spin of the gas, X g , within the halo is 
used instead, the discrepancy shrinks so that six of the 
eight galaxies have NFW Rd and i band results within 
10% of each othcr[3 As listed in Table [TJ the values 
of A are generally lower than X g for all halos except 

14 We note that the two that remain discrepant both show ev- 
idence for a downward "break" in the exponential light profile at 
large radii. 
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the two lowest mass galaxies. This trend for X g to be 
larger than A has been see n previously in simulations 
(jSharma fc Steinmetd [20051. and has been suggested as 
one reason why disks may have a smaller fraction of 
low angular mom entum tha n their dark matte r halos 
(jChen et al.ll2003D . Recently. iRoskar et all (|2010D used a 
similar simulation run with Gasoline to show that gas 
that enters the virial radius of a halo and cools toward 
the disk is torqued by the hot halo gas, so that the an- 
gular momentum of the disk gas becomes aligned with 
the hot halo. This is true even for gas that is initially 
counter-rotating with the disk, and hence is a powerful 
method to remove negative and low angular momentum 
material from the disk. The lack of low angular momen- 
tum material for the disk gas will lead to larger values 
of X g compared to A for the total halo. The disk stars, 
whose light we trace in the results of this paper, form 
from this cold gas with larger spin values, and thus X g 
will be a better predictor for the NFW Rd than A. 

The result that the disk baryons are lacking a low an- 
gular momentum component compared to the DM is not 
trivial. This result needs to be examined in detail. How- 
ever, this requires a full evaluation of the history of angu- 
lar momentum in these halos, which is beyond the scope 
of the current paper. We reserve such a study for future 
work. 

5. GROWTH OF THE STELLAR DISK SINCE Z = 1 

As discussed above, numerous observational studies 
have found that there has been a decline in the surface 
brightness of disk galaxies, by a magnitude or more, since 
z=l. Due to the fact that there is little evidence for a 
change in th e size function of disk galaxies over this time 
period (e.g., iKanwar et al . 1 120081 ) . the observations can 
only be interpreted as an upper limit to luminosity evo- 
lution. That is, if disks do not change in size, then all 
of the evolution must be in luminosity. However, due to 
the degeneracy in the magnitude - size plane, the option 
o f possible growth is l eft open. 

iBarden et al.l ()2005l ) investigated this disk growth. Af- 
ter showing that they, too, agreed with ~1 magnitude 
of surface brightness dimming back to z—1, they then 
used galaxy colors to derive stellar masses and considered 
the stellar mass - size relation for their galaxies. They 
found that this relation, unlike the magnitude - size re- 
lation, showed little or no evolution back to to z—1. Disk 
growth may still occur, but must occur in such a way that 
galaxies remain on the same stellar mass - size relation 
with time. They concluded that this was evidence for 
weak inside- out growth of galaxy disks, as did additional 
later studies (iTruiillo fc Pohlenl2005HTruiillo et aHl2006t 
iDutton et al ll2010bD . 

The right panel of Fig. [7] shows the evolution in the 
stellar mass - size plane for our simulated galaxies . 
In agreement with previous results (IBarden et al.l [20051 : 
iBrook et all 120061: iFirmani fc Avila-Reesd 12009ft . these 
galaxies are growing along a path that keeps them on 
the same stellar mass - size relation with redshift. 

The stellar mass shown in the right panel of Fig. [7] is 
for the disk only. To separate the disk growth from to- 
tal stellar growth of these galaxies (i.e., from the bulge 
and halo spheroidal stellar components) , a kinematic de- 
composition was done to identify disk stars at z=0. To 
identify disk particles, the galaxies are first aligned so 
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Fig. 9. — The surface mass density of the simulated disks as a 
function of redshift. The red points with error bars, connected by 
the solid black line, show the mean and standard deviation of the 
sample of nine at z = 0, 0.5, and 1. The mean varies by only 0.1 
dex between z=l and z=0. The dashed line shows the analytic 
prediction for the SIS model. The colored lines show the evolution 
of each individual disk galaxy. Colors are the same for individual 
galaxies as in Figs. [7] and [8] 

that the disk angular momentum vector lies along the 
2-axis. 3 z /3 C irc is calculated for each star particle in the 
galaxy, where J z is the angular momentum in the x-y 
plane, and 3 C irc is the momentum that a particle would 
have in a circular orbit with the same orbital energy. 
Disk stars are identified as those having near circular 
orbits, so that J z /J C i rc > 0.80 These z=0 disk stars 
were then searched for in their most massive progenitor 
at z=0.5 and z=l to find the mass of these stars that had 
formed at each step, with the results plotted in Fig. [Jj 

If galaxy disks are evolving approximately along the 
same stellar mass - size relation with time, then the 
surface mass densi ty of disks should s how little evolu- 
tion back to z—1 (IBarden et al.l [20051 : iSomerville et al.l 
120081) . That is, a galaxy with stellar mass of 10 10 M Q at 
z=l is roughly the same size (though maybe just slightly 
smaller) as a galaxy at z—0 with 10 10 M , yielding 
roughly the same surface densities. Fig. [H] tentatively 
confirms this result, showing that the surface mass den- 
sity of our population of disks evolves little with time, 
with surface mass density, £ defined as 

logYi — logM^^disk ~ ^loghs — log(2n). (5) 

The red points with error bars, connected by the solid 
black line, show the mean and standard deviation for 
the whole sample at z = 0, 0.5, and 1. The colored lines 
are the results for individual simulated disks with time. 
We note, however, that the dashed line shows the predic- 
tion for the su rface mass dens ity evolution for the sim- 
ple SIS model (|Mo et al.lll998D . The standard deviation 
within our small sample is large, so that the simulated 
galaxies are fully consistent with the predicted growth. 
While our sample is clearly too small to derive statisti- 
cal results, the change in the mean surface density of the 
population between z=l and z=0 is only 0.1 dex, half of 
the evolution predict ed by the SIS model. W e note that 
an updated model bv lSomerville et ail (|2008f) that incor- 
porates the evolution of NFW halos derived in N-Body 
simulations predicts only ^0.2 dex change in the stellar 
surface mass density back to z=l. This small evolution 

5 This criterion corresponds to an eccentricity < 0.2, which 
matches the eccentric ities observed in the Milky Way disk 
UNordstrom et al.ll200l 'l. 



Size-Magnitude Relation Evolution 



13 



in surface mass density is due to evolution along a stellar 
mass - size relation that changes little with time. 

5.1. The SFR - Mass Relation 

Due to the fact that the stellar mass - size relation 
shows little change with time, examining the evolution 
of galaxy disks at a fixed stellar mass is similar to ex- 
amining the evolution at a fixed size. However, because 
galaxies are growing with time, a galaxy that is, e.g., 3 
kpc in size and 10 M Q in stellar mass at 2=1 has in- 
creased in mass and size by z=0 and another galaxy has 
moved from a lower mass and smaller size into the same 
bin. Observations, and the results presented here, tell 
us that these two galaxies of same size at z=l and z—Q 
have a different luminosity: a galaxy at z=l with similar 
stellar mass and size to a z—Q galaxy will be brighter. 
Note that this is not equivalent to saying that an indi- 
vidual galaxy was brighter in the past. Instead, galaxies 
at fixed stellar mass were brighter in the past. Because 
they are the same stellar mass, the brighter 2=1 galaxy 
must then be forming more stars. 

Two of our galaxies can be compared in this way. 
Galaxy h603 has a stellar disk mass of 1.28xl0 10 M Q at 
z=0, while h258 has a stellar disk mass of 1.25xl0 10 
M at z=l. The SFR within the last lOOMyr for h258 
at 2=1 is 3.1 M /yr, while it is 0.6 M©/yr for h603 at 
2=0. The lower SFR for h603 is due to the fact that it 
resides in a less massive halo than h258 (see Tabled]). 
Even at z=l, h603 had a lower SFR than h258, of 2.1 
M Q /yr. This is because these galaxies follow a SFR - 
stellar mass relation (which also leads to their reproduc- 
tion of the observe d stellar mass - metallicity relation, 
Brooks et ai][2007l ). 

Fig. [10] shows the SFR - stellar mass and specific SFR 
(sSFR) - stellar mass relations for our simulated disks. 
The left panels of Fig. [TO] shows the sSFR (top) and SFR 
within the last 100 Myr (bottom) versus stellar disk 
mass. These left panels confirm the values listed in the 
paragraph above. The dwarf galaxies are included in this 
plot for full disclosure, but lie outside of the mass range 
observable at high z, unlike the higher mass galaxies. 
Two trends among the higher mass galaxies are note- 
worthy. First, the SFR of any individual galaxy tends 
to decrease slightly with time, but second, the entire 
population as a who le can be obse r ved to shift to lower 
SFRs and sSFRs l|Guo fc White! [2001 IDutton et all 
2010a; iBouche et all I2009Q . The decline in SFR and 
sSFR in has been observed extens i vely in populations 
of g a laxies dGavazzi fc Scodeggiol [ 1996; Boscll fet all 
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is one form of "downsizing" in galaxy evolution. 

The right panels of Fig. [10] also show SFR and sSFR 
for the simulated galaxy disks, but now as a function 
of total halo mass. This is included to demonstrate that 
although, e.g., h603 at z=Q and h258 at z=l have similar 
disk stellar masses, this does not correspond to similar 
virial masses for those redshifts (colors yellow and purple 
in the figures, respectively). The deeper potential well of 



the more massive galaxy (h258) leads to a higher SFR0 
In conclusion, attempting to compare galaxies that are 
a similar stellar mass at z=l and z=0 (and by the tight 
stellar mass - size relation, similar size) leads to inher- 
ently comparing galaxies of two different virial masses. 
This then leads to the observed difference in luminos- 
ity for galaxies of a similar size at z—1 and z—0, caused 
by the SFR - stellar mass relation, an d the correspond- 
ing stellar mass - h alo mass relation ()Guo et all 120101 : 
IMoster et aill2010bl ). 

5.2. Comparison to Observed Relations 

In Fig. [TTJ we reproduce some of the simulated scal- 
ing relations, this time with observed scaling relations 
directly overplotted for comparison. The left panel 
of Fig. [TTJ shows the low redshift relations for SDSS 
galax ies (solid line , with 2a scatter in dotted lines) 
from IDutton et all ( 2010b') . It is worth noting that 
IDutton et all ( 20101 ) performed bulge/disk decompo- 
sitions with the SDSS data, so that the sizes plotted 
are fits to purely exponential disk scale lengths (con- 
verte d from the ha l f light radius presented in their pa- 
per). IDutton et al.l (|2010b[ ) also use higher z data from 
the DEEP 2 survey to measure the change in size since 
z=0, finding ~0.1 dex of size evolution at a fixed mass 
back to z=l. In the mass range of overlap (M* > 10 9 
Mq), our simulated galaxies are entirely consistent with 
observed galaxies within the scatter. 

The central panel of Fig. [TTJ reproduces the stellar 
mass - SFR relati on, now with observational data from 
lElbaz et ail (|2007l) plotted for z=0 (solid line) and z=l 
(dashed line) , corrected from a Salpeter IMF to a Kroupa 
IMF. At both redshifts, the minimum mass observed is 
roughly 10 9 Mq. The right panel shows the same data, 
but now for the sSFR. The data point with error bars 
in the lower right corner of the righ t panel provides an 
estimate of the typical scatter in the lElbaz et al.l (|2007| ) 
data. While our galaxies are certainly consistent with the 
data within the scatter, and the general slope at each z 
is consistent with the observations, there is overall a ten- 
dency for the simulated galaxies to be more massive than 
the observations. We believe the discrepancy is in stellar 
mass rather than SFR, because these simulated galax- 
ies have both z=l and z=0 disk SFRs in good agree- 
ment with observed galaxies in the same halo m ass range 
(lConroy | fc Wechsle r 2009: Boissier et al1[2010l ). We note 
that the lElbaz et al.l (|2007l ) data is for total stellar mass, 
while our simulated results are for the disk only (bulge 
excluded). However, including the bulge mass will only 
exacerbate the discrepancy. For the massive galaxies, the 
discrepancy is probably due to the fact that these sim- 
ulations are known to for m too many stars (|Guo et all 
\2HWk . IMoster et al.l[2lTi0bh . 

We can learn something about when this overproduc- 
tion of stars occurs in the simulations by also examining 
the stellar mass - sSFR plot (right panel of Fig. ITT]) . 
While within the observed scatter, is appears that the 
sSFRs of the simulated disks are consistently lower than 

16 Note that although these two galaxies appear to have the same 
virial mass at the initial step, the properties of h258 have actually 
been measured at 2=1.25 due to the fact that it is undergoing a 
merger at 2=1. Galaxy h258 is 50% more massive than h603 at 

2 = 1. 
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Fig. 10. — The SFR and sSFR of simulated galaxies. Left panels: The bottom panel shows the SFR in the disk of each galaxy as a 
function of the stellar disk mass. The top panel is the specific SFR of the disk. Right panels: Same as the left panels, but now in terms of 
the total virial mass of the halo. The colors here are the same for individual galaxies as in previous figures. 



observational result s at both z=l and z=0 (lElbaz et al.l 
120071 INoeske et all I2007bb ISalim et al.l [2003). These 
plots demonstrate that these simulated galaxies must 
be forming too many stars at early times, at redshifts 
higher than plotted here (z > 1). Because the SFRs of 
the galaxies are in good agreement with observed galax- 
ies in the same halo mass range at both 2=1 and z=0 
(|Conrov fc Wechsleii |20Tj9|) . the fact that the sSFRs are 
lower than observed values requires that the stellar mass 
of these disks be too large by z=l. 

The story is slightly different for the two dwarf 
disk simulations. These t wo dwarf galaxi e s hay e 
been published previously in iGovernato et al.l (|2010f) . 
where it was demonstrated that their properties are 
a g ood match to observed dw arf disk galaxies (see 
also Sanchez- Janssen et al. 2010), including stellar mass 
(|van d en Bosc h et al.l 120011 ). However, the observed sS- 
FRs of low mass gal axies is on average higher than that 
of massive galaxies (|Salim et al.ll2007f ), indicating that 



low mass galaxies have been very inefficient in the past 
at turning gas into stars. Clearly, that is not the case 
with our two simulated dwarf disk galaxies, which have 
too low sSFRs compared to most observed dwarf galax- 
ies due to the fact that they form the bulk of their stars 
prior to z=l. Hence, while approximately the correct 
mass in stars is formed, the star formation is biased to- 
ward z > 1. We have verified that these two simulations, 
when run with stronger SN feedback, have burstier SF 
histories with their SF spread out more evenly down to 
the present time. This leads to approximately the same 
mass in stars formed by z=0, but with a more uniform 
SF history. This means that the stronger feedback dwarf 
runs have higher sSFRs than currently presented here, 
and sSFR values larger than the massive galaxy simu- 
lations shown here. Future work will include these new 
dwarf disks, which show similar stellar masses, rotational 
velocities, and sizes to those presented in this paper. 
Thus, in all of the simulated galaxies the star formation 
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Fig. 11. — Comparison of simulated disk galaxy scaling relations to observational results . Left: The stellar ma ss - size relation, reproduced 
from the right panel of Fig. [7] now showing (solid line) the z=0 observed relation of Dutton ct al. (2010b) plus the 2<r scatter in the 
observations (dott ed l ines). Their observed z=l relation would simpl y be shifted d ownward 0.1 dex in size. Center: The stellar mass — SFR 
relation, as in Fig. 1101 now showing the observed relations from Elbaz ct al. (2007) at 2=0 (solid line) and at z=l (dashed line). Right: As 
for the center pl ot, but now f or the sSFR. The x data point with error bars in the lower right corner of the panel represents roughly the 
2<t scatter in the lElbaz et all (120071 1 data. 



history is biased toward too much star formation at high 
z. By z=0, though, the dwarf galaxies have formed ap- 
proximately the observed amount of stars for their halo 
mass, while the more massive galaxies have overproduced 
the total amount of stars. It is suspected that this over- 
production of stars in the massive galaxy disks may be 
the result of the adopted low SF density threshold, which 
allows for most of the disk to be capable of SF at any 
given time, unlike the high density threshold which limits 
SF only to rare density peaks. However, as discussed in 
Section 2.1, it is currently computationally expensive to 
generate the more massive galaxies with the force resolu- 
tion required to resolve high density peaks where SF oc- 
curs, and no simulations of sufficiently high resolution are 
currently available to test this hypothesis. The adopted 
SF threshold is a compromise that allows stars to form at 
radii observed in real galaxy disks, and we have demon- 
strated that it leads to scale lengths in good agreement 
with observed galaxies. However, we note that even if 
our massive galaxies are forming too many stars, it is by 
roughly a factor of 3 at most, and at early times. This 
means that the more massive galaxies may shift up to 
0.5 dex to lower stellar masses in the stellar mass - size 
plane, but will still grow along the stellar mass - size re- 
lation, so that our conclusions are unaffected. Likewise, 
the B band magnitudes we present here are unlikely to be 
affected, as the bulk of overproduced stars must occur at 
z > 1. Since B band light is dominated by young stars, 
and our low z SFRs are in good agreement with observed 
galaxies in the same halo mass range, the magnitude - 
size relations presented in Figs. |4] and [5] are unlikely to 
change. Hence, the conclusions presented here about the 
origin of the evolution in the size - luminosity relation 
remain valid. 

6. CONCLUSIONS 

For the first time, a sample of very high resolution 
galaxy disks simulated within a fully cosmological con- 
text have been used to determine the evolution of the 
galaxy disk size - magnitude relation. Unlike previ- 
ous simulations, we determine disk scale lengths and 
magnitudes by fitting the stellar light profile rather 
than decomposing our simulations based on kinematics 



(and hence mass). We gene rate artificial surface bright- 
ness images using Sunrise (| Jonssonll2006t [Jonsson et all 
2010), in order to derive disk properties using methods 
used by observers. E xponential disk sc ale lengths were 
derived using Galfit (jPeng et al.ll2002D with redshifted, 
surface brightness dimmed i band images at z—1 and 0.5, 
and rest frame B band images at z=0. This allows for a 
direct comparison between simulated and observational 
results in the size - magnitude plane. 

Fig. @] demonstrates that these galaxies overcome past 
inabilities for cosmological disk simulations to retain an- 
gular momentum. Artificial loss of angular momentum 
is minimized due to high resolution and a physically mo- 
tivated SN feedback recipe, and hence we produce disks 
with realistic sizes compared to observed disks. Fig. [5] 
shows that the sizes are in agreement with observed disk 
sizes back to z=l. This success, and our previous suc- 
cesses in matching t he observed stellar m ass - metallicity 
relation with time (IBrooks et al.l 120071). a nd the Tully- 
Fisher relationship (jGovernato et al.H2009fl , indicate that 
these simulated disks reproduce the fundamental scaling 
relations for disk galaxies. Having established this, we 
can now use these simulated disks to investigate and in- 
terpret observational findings. 

We use the simulated disk properties to investi- 
gate the evolution in the magnitude - size relation 
with time. Observations have generally concluded 
that there has been 1 to 1.5 magnitudes of sur- 
face brightn e ss dimming since z =l dSchade et al.lll99€ 
Roche et all 119981: ILillv et all 119981: iBouwens fc Silt 



2002; Ravindranath ct al. 2004; Truiillo & Pohlcn 2005; 
Bardenet al. 2005; Melbourne et al.ir2007; Kanwar et alj 



20081) . We find similar evolution for our galaxies (see 
Fig. Combined with a lack of evidence for an 

evolving size functio n for galaxy disks (ILillv et al.lll998h 
iRavindranath et al.l l2004t iKanwar et alj 120081 ). and lit- 
tle to no growth i n galaxy disks sizes bac k to z—\ 
(|Barden et all [2001 iTruiillo fc PoTIenl [20051 ). the ob- 
served surface brightness evolution can only be inter- 
preted as an upper limit to luminosity evolution if there 
has been no change in galaxy disk sizes. 

We have shown that we 1) are in good agreement with 
the observed magnitude - size relation with time and 
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the evolution in surface brightness, but 2) a number of 
the simulated disks are clearly undergoing a large change 
in size (Figs. [7] and [8]). Our ability to follow individual 
galaxies with time, unlike the observations, allows us to 
interpret the evolution in the magnitude - size plane. We 
find that the evolution is dependent on mass, with our 
massive galaxies undergoing larger changes in size than 
magnitude. 

Our dwarf galaxies undergo the least change in 
size, though there is no immediate theoretical expec- 
tation that this shou ld be the case. Recent results 
(|Fakhouri et al.l [20101 ) show that the rate of growth of 
halos has little dependence on mass, and thus low mass 
galaxies should grow just as much in size as their more 
massive counterparts. We conclude that a larger sam- 
ple is needed to study if low mass galaxies truly grow 
less than more massive galaxies. On the other hand, 
the evolution in magnitude for dwarf disk galaxies is de- 
pendent on their bursty SF history, but an individual 
dwarf can undergo significant dimming of at least 2 mag- 
nitudes. A dichotomy of evolution with mass has been 
seen in observations, with low mass galaxies exhibiting 
significantly more dimming since z=l than more mas- 
sive g alaxies (e.g.. [Melbour ne et al.ll2007t iKanwar et al.l 
f2008ff . 

The halo properties for these galaxies (A, c, md,jd) pre- 
dict a scale length that is generally shorter than the scale 
lengths derived from the simulated light profiles. Fig. |3] 
demonstrates that the light and mass distributions yield 
scale lengths in good agreement, so the discrepancy be- 
tween Rd predicted for NFW halos and from our light 
profiles cannot be due to a discrepancy between light 
and mass profiles. Rather, using A s , the spin of the 
gas in the halo, brings the predicted and measured scale 
lengths into decent agreement. Thus, the angular mo- 
mentum distribution of the gas from which the stars are 
forming has been modified from that of the dark matte r 
(see also iSharma fc Steinmetzl l2005t IChen et al.l [2003). 
iRoskar et al.l (|2010f ) have demonstrated that accreted gas 
is torqued by the hot halo after entering the virial radius, 
potentially preventing counter-rotating and low angular 
momentum gas from reaching the disk. A full examina- 
tion of how this process occurs is beyond the scope of 
this paper and left for future work. 

Having determined that these galaxies are growing in 
size with time, we can investigate this gr owth in stel- 
lar m a ss. Like previous t heore ti cal mode ls (Broo k et al. 
20061: iSomerville et al.l 120081 : iFirmani fc Avila-Reese 
2009; IDutton et aUl2010blT we find that these galaxies 
grow in such a way as to stay along approximately the 
same stellar mass - size relation with time. Hence, at a 
fixed stellar mass, there is little change back to z=l in 
the size of galaxies. 

Due to the fact that our larger galaxies have grown 
since z=l, we do not have the ability to directly com- 
pare results at a fixed size. In particular, we do not have 
a population of simulated disks at z=l with scale length 
greater than 3 kpc, though these larger disks do exist 



observationally. Presumably we might generate these 
larger disks if we simulate higher mass halos. How- 
ever, it is not clear from this study how these larger 
disks would evolve to z—0. Observationally, the num- 
ber density o f disk galaxies shows no evidence for change 
back to z=l ([Lilly ■ et al.lll998tlRayindranath et al.ll200l 
iSargent et al.ll2007tlKanwar et al.ll2008|) . In coniunction 
with the apparent lack of evolution in the size function 
of disks, this fact has been interpreted to mean that ei- 
ther galaxies are done growing by z=l, or that the rate 
of destruction of galaxies at given size must be equally 
matched by the rate of growth of galaxies into that size. 
The fact that our galaxies are growing, while still match- 
ing the size-magnitude evolution, suggests that galaxies 
are moving into and out of a given bin in size during this 
time interval. It is beyond the scope of this work to an- 
swer what happens to these larger galaxies, though we 
may speculate. There are two possibilities that arise due 
to mergers; perhaps the larg er disks undergo mergers and 
become early type galaxies ()Bell et al.ll2007D . or perhaps 
mergers cause them to become bulge dominated disks. 
As the majority of size-magnitude evolution studies have 
focused on pure disks or disks with small bulges, this 
would remove them from the sample being considered 
(jSargent et al.ll2007ft . 

Finally, we demonstrate that the simulated galaxy 
disks follow SFR - mass and sSFR - mass relations 
similar to observations. While the growth of dark 
matter halos is nearly self-similar with mass, observa- 
tions demonstrate that some proc ess must break this 
self-similarity in the baryons (e.g.. iBenson et al.l 120031 : 
IConrov fc Wechslerl [2001 ISchave et al.l 12010ft . In the 
mass range of galaxy halos presented here, SN feedback is 
the process that regulates SF as a function of halo mass, 
and likely drives gas outflows that vary as a function 
of mass (though outflow rates are left for future work). 
The regulation of SF as a function of mass leads to the 
reproduction of the observed st ellar mass - metallicity 
relationship (|Brooks et al.l 12007). and the SFR - mass 
relation shown here that is responsible for the reproduc- 
tion of the size - luminosity relation back to z=l. A 
galaxy at a fixed M* at z=l will have a higher SFR than 
a z=0 counterpart, and thus a higher luminosity as well. 
In conjunction with the weak evolution of the stellar mass 
- size relation, this tells us that the difference in lumi- 
nosity observed at a fixed size between z=l and z=0 is 
due to the fact that the z=l galaxy has a higher SFR 
and luminosity. 
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